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We study static, spherically symmetric black holes in the recently proposed mas- 
sive gravity theory with two dynamical metrics that is maintained to be ghost-free. 
These solutions possess a regular event horizon common for both metrics, the val- 
ues of the surface gravity and Hawking temperature calculated with respect to each 
metric being the same. The ratio of the event horizon radii measured by the two 
_ . metrics is a free parameter that labels the solutions. We present a numerical ev- 

(— 1 1 idence for their existence and find that they comprise several classes. Black holes 

within each class approach the same AdS-type asymptotic at infinity but differ from 
each other in the event horizon vicinity where the short-range massive modes reside. 
In addition, there are solutions showing a curvature singularity at a finite proper 
distance from the horizon. For some special solutions the graviton mass may become 
effectively imaginary, causing oscillations around the flat metric at infinity. The only 
asymptotically flat black hole we flnd - the Schwarzschild solution obtained by iden- 
tifying the two metrics - seems to be exceptional, since changing even slightly its 
horizon boundary conditions completely changes the asymptotic behavior at infinity. 
We also construct globally regular solutions describing 'lumps of pure gravity' which 
show the same far field behavior as the black holes. However, adding a matter source 
gives asymptotically fiat solutions exhibiting the Vainstein mechanism of recovery of 
General Relativity in a finite region. 
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I. INTRODUCTION 



A possibility that gravitons may have a mass |l| was thought for a long time to have 
mainly an academic interest. However, the recent observational evidence 2! suggests taking 
this idea more seriously, because it could provide an explanation for the current acceleration 
of our universe jsl. This has produced an increase of interest towards massive gravity theories 
(see for a recent review). 

Such theories are known to have serious theoretical difficulties - the absence of a smooth 
massless limit ^, presence of the ghost state in the spectrum jol, and the very low ultra- 
violet cutoff [7|. However, it seems that remedies may exist for some or perhaps for all 
of these problems. For example, the van-Dam- Veltman-Zakharov (vDVZ) discontinuity in 
the massless hm.t seems to be eurable by the Vaiusteh. ,„eeha„,s,„ g as was recently 
confirmed by the explicit calculation [9|. In addition, it seems that the presence of the ghost 
and the absence of uniqueness may cure each other, since among many possible theories of 
massive gravity there could be one that is ghost-free. 

A special massive gravity model which is the only one that is ghost-free in the decoupling 
limit was recently discovered by de Rham, Gabadadze, and ToUey (RGT) [l^. A number of 
reports have then maintained that this theory should be completely free of the ghost 



although the opposite claim has been advocated as well [12|. Setting aside this apparently 
difficult issue, one should say that the theory of jlO|, called below RGT model, is certainly 
distinguished. For example, some of its equations of motion 'miraculously' become algebraic, 
instead of bein g di fferential. This model has received a lot of attention, in particular its static 



solutions 



13|, IJ] and time-dependent cosmological solutions 16| have been studied 



The RGT model is a bimetric theory containing a dynamical metric gfj_u and a non- 
dynamical flat reference metric /^^. Quite recently, its bigravity generalization was proposed 
by Hassan and Rosen (HR), within which the metric /^,^ is also dynamical, so that the theory 



contains two gravitons - one massless and another one massive [17|. Moreover, it has been 
maintained that the generalized model remains ghost-free jl8|. This theory contains the RGT 
model as a special limit, but it has richer dynamics, since it contains more degrees of freedom. 
For example, it reproduces General Relativity when the two metrics are identifled, whereas 
otherwise the massive degrees of freedom are present. In this respect the theory resembles 
the standard gauge fleld theories, where the massive states are generated by assuming a 
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non-zero value for the Higgs field, but one can switch off the mass by setting the Higgs field 



to zero. Cosmological solutions in the HR theory were studied in 
solutions were obtained in |2l| . 



19[, 20| while some static 



In the present paper we carr y o ut a detailed analysis of static, spherically symmetric 



black holes in the HR theory of 



17| . Our motivation is provided by the fact that the only 



known up to now black holes in massive gravity are all described by the Schwarzschild- 
(anti)de Sitter type metrics. These are exact solutions for which the metrics g^i, and f^^, are 
not simultaneously diagonal, they were discovered long ago in a particular massive gravity 
model 



, but later similar solutions were found in many other massive gravity theories 
15|. It is, however, unclear if other, asymptotically fiat, say, black holes can exist 



in massive gravity. 

The two metrics can also be chosen to be simultaneously diagonal, but this leads to rather 
complicated equations. These equations admit asymptotically fiat (numerical) solutions 
when a regular matter source is included |9|, but it is unclear if their black hole solutions 
exist as well. The recent results seem to exclude such a possibility at least in the bimetric 
models where one of the metrics is non-dynamical and fiat. Specifically, it turns out that 
when one numerically integrates the field equations starting from infinity towards an inner 



region, one does not find a regular event horizon but a naked singularity instead 23[. In 
addition, a geometric analysis of 2^ shows that if a regular horizon is present, then it should 
be common for both metrics, which is impossible if one of them is flat. 

At the same time, asymptotically flat black holes certainly exist in the bigravity theory, 
where both metrics are dynamical. These are the usual vacuum black holes obtained by 
setting g^i, = f^^. It is then natural to try to find their generalizations to the case where 
g^u 7^ ffjLu and the massive degrees of freedom are excited. Such solutions could describe 
black holes surrounded by a cloud of massive hair. We therefore study in what follows black 
holes in the theory of |l7|, focusing in the main part of the paper on the case where both 
metrics are diagonal, the opposite case being considered in the Appendix. 

We assume a regular even horizon and starting from it integrate numerically the field 
equations towards infinity. We find that the horizon is common for both metrics and that 
the values of the surface gravity and Hawking temperature calculated with respect to each 
metric are the same, which agrees with the conclusions of Ref. 2J|. The ratio of the event 



horizon radii measured by the two metrics is a free parameter that labels the solutions, 
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so it is sufficient to vary only this parameter to construct all possible black holes. We 
then proceed to classify the solutions and find the following generic types of behavior: the 
solutions approach at infinity either the anti-de Sitter (AdS) metric, or special metrics that 
we call U,a backgrounds (see Eq.f l4.7p below), or they develop a curvature singularity at a 
finite proper distance from the horizon. There are many solutions which run to the same 
AdS or f/, a asymptotic at infinity but differ from each other in the event horizon vicinity, 
where a short massive 'hair' is localized. On the other hand, the asymptotically flat solution 
of Schwarzschild seems to be exceptional, since changing even slightly its horizon boundary 
conditions immediately changes the asymptotic behavior at infinity. As a result, we do find 
black holes with massive hair, but they are not asymptotically flat. It seems that asymptotic 
flatness is incompatible with the presence of a black hole horizon when the massive degrees 
of freedom are excited. 



For the sake of completeness, we also briefly study solutions without a horizon which 
describe 'lumps of pure gravity'. They have a regular center, but they are also not asymp- 
totically flat and show in the far field the same behavior as the black holes. It seems they 
can be obtained from the latter in the limit where the horizon shrinks to zero. However, 
adding a matter source we do find asymptotically flat solutions which exhibit the Vainstein 
mechanism of recovery of General Relativity at finite distances. This suggests that the 
mechanism needs a matter source and so does not work for pure vacuum systems like black 
holes. 

In the following three sections we describe the bigravity theory of its equations 
of motion, the reduction to the spherically symmetric sector, and simple exact solutions 
that we call 'background black holes' and 't/, a backgrounds'. Sec.V explains our procedure 
of numerical integration of the equations starting from the horizon. In Sec. VI we study 
solutions which are parametrically close to the analytically known background black holes. 
The generic parameter values are considered in Sec. VII, with some limit cases analyzed in 
Sec. VIII. In Sec. IX we present globally regular solutions without a horizon, while Sec.X. 
contains concluding remarks. Finally, the Appendix describes exact black hole solutions 
with non-simultaneously diagonal metrics. 
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II. THE BIGRAVITY THEORY 



The theory of [17| is defined on a four-dimensional spacetime manifold Ai spanned by 
coordinates x'^ and equipped with two metrics g^y{x) and f^p{x). Their kinetic terms are 
chosen to be of the standard Einstein-Hilbert form with the couplings denoted by G and 
tan^ 7] G. We parameterize the action as 

S=--^ I f^i^ + m^cos^r/A^t) v^^'x- \ [ U^/^ d^x + S^^) , (2.1) 
SttG J \2 J lovrtan r/G j 

where R and 7^ are the Ricci scalars for g^^y and /^,^, respectively, while S-^a describes ordinary 
matter (for example perfect fiuid) which is supposed to directly interact only with g^^y. The 
interaction between the two metrics is defined by 

Ant = ^(^' - K^^) + I ^,u,.e-^'''K^^K^f,K^, + I e,y,^e-f'^<'K^^KmKl (2.2) 
with = 6^ — 7^j^ and 7""^ defined by the relation 

r.i''. = g'^'U (2.3) 

where g'^'^ is the inverse of g^^. Apart from the gravitational coupling G the theory contains 
three parameters 77, C3, C4 as well as the graviton mass m. In the limit where 77 — )■ and f^^, is 



fiat the theory reduces to t 



le RGT model of 



10] . The above parametrization of the theory 



was also considered in Ref.|l9| (with tan 77 — )■ 7], mcost] — )■ m, Kj^ — —Kj^, C3 — )■ — C3). 
In order to calculate 7^^ defined by (12. 3p . it is convenient to introduce two tetrads 
and defined by the conditions g^'^ = rj^^e'^e'^ and = rjAB^^^y with rjAB = 

diag(l, —1, —1, —1). These tetrads are defined up to the local SL{1, 3) x SL{1, 3) rotations, 

which freedom can be used to impose the conditions 

e^^ws^ = e^WA/x (2.4) 
with ua^j. = Vab^I^j which insures that 

7^ = e^;u;^ (2.5) 

Varying the action gives the field equations (see jl9| for details) 

Gi = cos' r/ + MGT^'^^i , (2.6) 
g{ = m'siYi\Tl, (2.7) 
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where and are the Einstein tensors for g^j^^, and /^,^, respectively, while 



Tj: = rf-5,^£„,, ^' = -73f<' (2.8) 



with 



< = (7: - 3)7^ - 7^A + I eA,..e"''^'^7^i^;i^; + | 6,,..6"/'^^^i^^^i^;i^,- . (2.9) 
The Bianchi identities for the left-hand side of Eg. (12.61) imply the conservation condition 

Vpr( = o, (2.10) 

(9) 

where Vp is the covariant derivative with respect to gfj,,y. In view of the diffeomorphism- 

{/) 

invariance of the interaction term S'int = / jC.mty/—g d^x, the conditions Vp = then 
follows identically. The matter energy-momentum tensor is conserved independently, in 
view of the diffeomorphism-invariance of the matter action S'(m), 

v^r(-)^ = o. (2.11) 

If the matter source vanishes, then choosing the two metrics to be the same reduces the field 
equations to those of vacuum General Relativity, 

9,u = f,u, T^'''\ = G^ = ^,^ = 0, (2.12) 

since one has in this case = = 0. 

In what follows we shall be considering solutions of equations (12. 6p . (12. 7p within the 
spherically symmetric sector. Introducing the spherical coordinates = (t, r, t?, y?), the 
most general expressions for the two tetrads are [19 1 

_ I d _ d _ 9 

""""gat' 9^' ""^'rW ""''R^d^' 

= adt + cdr, = -cQN dt + bdr, oj^ = Ud-d, oj^ = UsmMip, (2.13) 

where Q,N,R, a, b, c, U are functions of r. The corresponding metrics read 

dv'^ 

g^^dx^'dx" = Q^de -j^- R\d^^ + sin^ Mip^) (2.14) 

and 

f^^dx^'dx" = (a^ - c^Q^N^) dt^ + 2c(a + hQN) dtdr - (6^ - c^) dr"^ - U'^{d§'^ + sin^ M^"^), 

(2.15) 
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while 

/ a/Q c/Q Q \ 
-cQN'^ 6iV 
U/R 
\ U/R J 
It is now straightforward to compute C\n f an d the tensor r/f defined by (12. 9p (the exphcit 



7^ 



(2.16) 



expressions are given in the Appendix in 



19|). 



It is still possible to reparameterize the radial coordinate to impose the gauge condition 
R{r) = r so that the metric g^^ becomes 



g^udz'^dx" = Q'^dt^ 



Ar2 



-r'(cir + sin' 7? V)- 



(2.17) 



Since its Einstein tensor is diagonal, so should be the energy-momentum tensor. Therefore, 
one has to have T° = 0, which requires that r° = 0. Using (12. 9 p one finds 



(r {2U - 3r) + cg (3r -U){r-U)+ C4 (r - Uf) 



f2.18) 



and for this to vanish, one should either have c = 0, or choose c 7^ but set to zero the 
expression between the parenthesis. The latter case is studied in the Appendix, where it is 
shown that the metric g^y will be described in this situation by the Schwarzschild-(anti)de 
Sitter solution. More general metrics are obtained by choosing c = 0. 



III. FIELD EQUATIONS 

If c = then the metric f^y in (I2.15P becomes diagonal. Assuming that all its coefficient 
depend only on r and introducing a new function Y via h = U'/Y with ' = d/dr the metric 
reads 

TT/2 

f^.^dx^'dx" = a^de -y^dr"^- U'^{d^'^ + sin^ M^"^). (3.1) 

We thus have two static and spherically symmetric metrics (I2.17P and (13. ip which contain 
5 functions of r: Q, A^, a, Y , U. The Einstein equations (r2.6|] . (l2.7p reduce then to the 
following system, 

G° = cos' r/T° + p, 
G;: = m' cos' r/T;-P, 

r>(^ 2-2 -T-O 

= m' sin' r/r;, (3.2) 
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(9) „ 

which should be supplemented by the conservation condition Vp = 0, 

(t;)' + ^ (t; - T") + ^{T^ - t;) = o. (3.3) 



(/) 

As was said above, it is not necessary to require in addition that Vp = 0, since this 



condition follows from fl3.2p .f l3l3l) . It is assumed in the above equations that the matter 
source is of the perfect fluid type, 

SvrGT^"'^^ = diag[p(r), -P(r), -P(r), -P(r)] (3.4) 

whose conservation requires that 

P' = -|(p + P). (3.5) 

The 5 equations f l3.2p .f l3l3l) explicitly read 

2NN' . 2 f N , \ , , 

H 5 h m cos 7] [ ai — U + a2 + p = 0, (3.6a) 



2N^Q' - 1 



, ^2 , m cos 7/ (^ai - + aaj - P = 0, (3.6b) 

{F^ - 1 + sin^ r] a^jNU' + 2UNYY' + sin^ t] Ya^ = 0, (3.6c) 

{a(F2 - 1) + sin^ r/ a^jU' + 2t/r2a' = 0, (3.6d) 

agf/' + «7a.' = 0, (3.6e) 

where the following abbreviations have been introduced, 

2(C4 + 2C3-I)f/ (c4 + C3)t/2 



«! = 3 - 3C3 - C4 + 

2(3 - C4 - 3C3)f/ (C4 + 2C3 - 1)U^ 



«2 = 4C3 + C4 - 6 



"3 = C4f/^ - 2(C3 + C^yU + (C4 + 2C3 - 1^ , 

a4 = (3 - C4 - 3c3)r2 - + c^)U^ + (4c3 + 2c4 - 2)rf/, 

"5 = [(a - Q)c4 - Qc^V + [2(2Q - a)c3 + (Q - a)c4 - g]r[/, 

+ [(2a - 3g)c3 + (a - g)c4 + 3g - ay, 
ae = Q'N[{3c3 + C4 - Sy + (2(1 - C4 - 2c3))f/r + (04 + 03)^7'], 

+ 2g(r - A^)[(3 - C4 - 3c3)r + (04 + 2c3 - l)U], 

+ 2a(iV - F)[(l - C4 - 2c3)r + (04 + C3)f/], 
ay = F[(3 - C4 - 3c3)r2 + 2(c4 + 2c3 - l)Ur - (04 + C3)U^]. (3.7) 



9 



Eliminating the term Q' in ag with the use of (13.6bp . the coefficients ai, . . . , ay depend only 
on Q, iV, a, y, U but not on their derivatives. 

The equations admit the scale symmetry that maps solutions to solutions, 

N{r)-^N{\r), F(r) ^ F(Ar), U{r) ^\u{\r), Q{r)^Q{\r), 

A 

a{r) -> a(Ar), m — — , rj ^ t], C3 — C3, C4 — C4. (3.8) 

A 

In the following few sections, until Section IX, we shall set p = P = 0. 



IV. SIMPLEST SOLUTIONS 

Some exact solutions of equations (I3.6ap - fl3.6ep can be obtained. 

A. Background black holes 

Let us choose the two metrics to be conformally related, 

Uu = C^g,., (4.1) 

with constant C. Equations (I3.6ap - fl3.6ep will be fulfilled if C satisfies the algebraic equation 

(C - l)P(C) = 0, (4.2) 

with 

P{C) = (c3 + C4)C2 + (3-5c3 + (e-2)c4)C+(4-30c3 + (l-20c4-6 + ^^^^^±^i^, (4.3) 

(here ^ = tan^?]), while g^y is the Schwarzschild-(anti)de Sitter metric, so that 

Y^ = N^ = l- — -^r\ U = Cr, a = CQ, Q = qN, (4.4) 
r 3 

where q is an arbitrary constant related to the time scaling symmetry, while 

A(C) = m2cosS/(l-C){(c3 + C4)C2 + (3-5c3-2c4)C + 4c3 + C4-6}. (4.5) 

One root of Eq.f l4.2p is C = 1, in which case A = and we obtain the Schwarzschild solution. 
Depending on values of the parameters C3, C4, 77, equation (14. 2 p can have up to three more 
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real roots, for which A(C) does not generically vanish. For example, for m = 0.1, 77 = 1, 
C3 = 0.1, C4 = 0.3 Eq. fi4.2p has altogether four real roots, C = {C^}, 

{Ci, C2, Cs, C4} = {1; -0.6458 ; 2.6333 ; -8.5566}, 

^^^^ = {0; -3.0559; -1.1812; +21.5625}. (4.6) 

This gives the Schwarzschild (S) solution for C = Ci, Schwarzschild-anti-de Sitter (SAdS) 
solutions for C = C2, C3, and the Schwarzschild-de Sitter (SdS) solution for C = C4. We 
shall call these solutions background black holes, since below they will be considered as 
reference backgrounds - the starting point for studying more general solutions. 

It is worth noting that these solutions are not the same as those in the Appendix, where 
Qfj^iy is also SdS or SAdS but g^^, and f^^, are not proportional. 



B. [/, a backgrounds 

Another class of solutions can be obtained by setting U, a to constant values: 

(2Af 4 1 

(1 - 2c3 - C4)f/^ + (3c3 + C4 - 3)Ur + (2 - - C3 - - c^y 
r 3 3 

Q rn? cos^ rj T dr w? sit? t] r dr 

N^"" 2 X, ^iV^ ' 2U Jr, N ' 

where = (C4 — 3 + 3c3)x^ + 2(1 — 2c3 — C4)Ux + (03 + 04)?/^ and M, ri, r2 are integration 
constants. For r — > 00 one has ~ N'^ ~ y ~ and the metric g^i, approaches in 
the leading order the (anti)de Sitter metric, but the subleading terms are different. The 
metric ff^^, is actually degenerate, since frr = W^/Y"^ = 0. However, such solutions will 
describe below the asymptotic behavior of other, more general solutions for which U, a 
become constant only for r — t- 00 so that /„ vanishes only asymptotically. The proper 
distance till infinity {U' /Y)dr and the proper volume of the 3-space are then finite, so 
that with respect to /^^ the spacetime is spontaneously compactified. 



V. BOUNDARY CONDITIONS AT THE HORIZON 



Since we could not find other solutions of equations f l3.6ap -( !3.6ep in a closed analytic 
form, we wish to integrate these equations numerically. As a first step, one should resolve 
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the equations with respect to the derivatives, and our procedure is as follows. First of all, 
taking the ratio of Eqs. fl3.6dp and fl3.6ep leads to the algebraic relation 

a{Y^ - 1) + sin^ T] a5 2UY^ 



which can be expressed in the form 



Using this, Eqs. (l3.6ap . (13.6bp . (I3.6cp and ( I3.6ep reduce to 



where J-^ = J-fc(r, TV, F, ?7, m, r/, C3, C4). Eqs. (l5.2p . (I5.3cp . (15.3dp together imply that 



(5.1) 



^ = ^(r,iV,F,f/,m,r/,C3,C4). (5.2) 



N' = J^^U' + J^2, (5.3a) 

Y' = J^^U' + :Fi, (5.3b) 

g' = -Fs Q , (5.3c) 

a' = TQQU', (5.3d) 



J^' = TqU' - 77^ , (5.4) 
and calculating the derivative on the left gives 

drT + d^TiTx U' + T2) + 9yJ^( U' + Ti) + dijTU' = T^U' - TT^, (5.5) 
which can be resolved with respect to [/', 

T% — T\d^T — J^-idyJ^ — duJ^ 
Injecting this into (I5.3ap and ( ]5.3bp finally yields a system of 3 coupled equations, 

N' = VN{r, iV,F,t/,m,r/,C3,C4), 
Y' = VY{r,N,Y,U,m,r],cs,c,), 

U' = VU{r, N, Y, U, m, r/, C3, C4), (5.7) 

where the explicit expressions of the functions on the right are somewhat lengthy, so that 
we do not write them down explicitly. When a solution of these equations is found, then Q 
is obtained by integrating equation (I5.3cp . and finally a is obtained from (15. 2p . 
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We wish to study black hole solutions of Eqs. (15. 7p and so assume that there is an event 
horizon at r = r/j where Q{rh) = N{rh) = 0. For the horizon to be non-singular and 
non-degenerate, and A^^ should both have simple zeros at r = Vh.. Next, the analysis of 
Eqs. (15.71 ) shows that if A^^ has a simple zero at r = r/^ then Y"^ should have a simple zero 
at this point too. At the same time, U can assume at the horizon any finite value. We 
therefore assume the power-series expansions near r = r^, 

N' = J2 «n(^ - rhT, y' = Yl ^"(^ - ^'^)" U = urf, + J2 Cn(r - r,)". (5.8) 

n>l n>l n>l 

Injecting this to (15. 7p . all coefficients a„, 6„, c„ can be expressed in terms of u, ai where u is 
arbitrary while ai satisfies A a\+B ai+C = 0, where A, B, C are functions of u, Vh, m, rj, Ci, C2. 
There are two solutions for Oi, 

1 



ai = af{u) = —{-B± VB^ - AAC), (5.9) 

which give rise to two different local solutions (15. 8p with different values of {N"^)', (^^)', U' 
at r = Th- Using these solutions in Eqs. (15.21) . (15. 3cp gives 

n>2 n>l 

where c„ and dn are expressed in terms of ai and u, while g is an integration constant that 
reflects the possibility to rescale the time coordinate for both metric simultaneously. 

At this point it is interesting to compare our results with the geometric analysis of 



bimetric theories of Ref. 



24| . Let us consider C, = d/dt, the timelike Killing vector for the 



both metrics. Its norms calculated with respect to each metric, (^, C,)g = and f — ^ i 
both vanish at r = r/^ so that the horizon is the Killing horizon for each of the two metrics 
at the same time. The surface gravity for each metric is the horizon value of 

= -I g'^-g., V, C V„ e = lim Q'N" = ] q'a, , 

4-4 't e e - m„ ) 1 |^ . (5.n) 

The Hawking temperature for each metric is obtained in the standard way by passing to the 
imaginary time and requiring the absence of conical singularity. This gives 
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Now, using the explicit expressions for the coefficients ai, bi, Ci, di (we do not write them 
down in view of their complexity) the ratio of the two surface gravities evaluates to one, 



K 



dibi 



(5.13) 



for any u and for both signs in (15. 9p . Therefore, the surface gravities and the Hawking 
temperatures are the same. All this agrees with the conclusions of Ref.j2^ that the two 
Killing horizons and their surface gravities should coincide. 

Returning to our analysis, the black holes solutions are obtained by numerically extending 
the local solutions (15. 8p . (I5.10p towards large r. It follows that the black holes are determined 
by the value of m = U{rh)/rh - the ratio of the event horizon radius measured by /^j, to 
that measured by g^^,^ - as well as by the choice of sign in (15. 9p . Depending on the latter, 
we shall say that the solutions belong either to the upper or to the lower branch. Under the 
scale transformations (13. 8p the ratio u = U{rh)/rh stays invariant while — ?■ t^/X^ and this 
can be used to set r/^ = 1, which condition will be assumed from now on. 

It is then relatively easy to explore the structure of the solution space, since it is enough to 
integrate the equations for various values of m, separately for each branch. The integration 
constant q in (I5.10p is not an essential parameter and can be fixed afterwards, when the 
global solutions are already known, for example by requiring that Q'^ /N'^ — 1 as r -> oo. 
Let us choose some values of the theory parameters, for example the same as in Eq. (l4.6p : 
m = 0.1, ?7 = 1, C3 = 0.1, C4 = 0.3. There is nothing special about these values, since varying 
them changes the solutions smoothly, without changing their qualitative structure. We then 
use the numerical routines of j25| to integrate the equations starting from r = r/, = 1 towards 
large r. 

To begin with, for m = C with C = {Ck} given by Eq. (14.60 . the numerical solutions 
should reproduce the background black holes, for which U'{r) = C. Indeed, for the four 
values of C in Eq. (l4.6p we find numerical solutions with U'{r) = C and N,Q,Y coinciding 
with those in Eq. (l4.4|) . The solutions for C = {Ci,C2} belong to the upper branch, while 
those for C = {C3, C4} belong to the lower branch. 
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VI. HAIRY BLACK HOLES 

As a next step, we choose u = C + 5u where 5u is smalL It is then natural to expect 
the solution to be a slightly deformed background black hole, the deformations being due 
to the 'massive hair' present for generic values of u. In what follows we analyze these 
expectations and find that the deformations are indeed small, but in general only within a 
finite neighborhood of the event horizon and not necessarily for all values of r. 

A. Deformations of the asymptotically fiat black hole. 

The only asymptotically flat solution among the background black holes is the Schwarzschild 
solution fl4.4p with C = Ci = 1. We consider its deformed version for n = 1 + 5u. It turns 
out that choosing 5u to be negative does not give anything, since the local solutions ( 15. 8p 
become then complex-valued. However, for 5u positive and small enough, for example 
5u = 10~^, the deformed solution exists and stays very close to the Schwarzschild solution in 
a large vicinity of the horizon, for r < rinax(^)- However, for larger values of r it completely 
changes its structure, since the Q, N, Y amplitudes then grow rapidly, while a, U approach 
finite asymptotic values (see Fig.l), so that the whole configuration approaches one of the 
U, a backgrounds (14. 7p . 
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Figure 1. Left: solution profiles for u = 1.01. Right: U'{r) for several values of u. Here and in 
Figs. 2, 3, 4 below one has m = 0.1, r] = 1, cs = 0.1, C4 = 0.3. 

It is instructive to consider the function U'{r), which is equal to one everywhere for u = 1. 
For u > 1 it stays very close to one for r < r^i^{u) but at r ~ r^i^{u) it suddenly drops down 
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and after a couple of oscillations tends to zero at infinity. The value r^g^^{u) increases when 
u decreases, and in the limit n — )■ 1 one has rjnax(^) — )■ oo so that the solution approaches 
the Schwarzschild metric for any finite r. However, since the boundary conditions at infinity 
for u > 1 are not the same as for u = 1, the convergence in the limit m — 1 is only pointwise 
and not uniform. 

The conclusion is that exciting the massive degrees of freedom around the Schwarzschild 
black hole produces deformations which stay small close to the horizon but inevitably become 
large at infinity, thus destroying the asymptotic flatness. 

The above conclusion is supported by the following counting argument. Let us require 
the solution to be asymptotically fiat. Then for r — j- oo one should have = 1 + 6N, 
Y = 1 + SY, U = r + 6U where the variations are small. Inserting this into Eqs.( l5.7() and 
linearizing with respect to 6N, 6Y, 6U and keeping only vanishing at infinity modes gives 



N = 1 + B cos^ ?7 e ™^ U = r + B — e "'"^ 

A siv? ri ^ ^ l + mr 
Y = l is sm rj e , (6.1) 

where A, B are integration constants. The other two metric amplitudes then read 

g = l-^i^ + ^^^e— a = l-^i^-'^^^"'\— . (6.2) 

ry ry ry 

This asymptotic solution is the superposition of the long-range massless Newtonian mode 
and the short-range massive VdVZ mode jsf. 

Suppose that one wants to find black hole solutions with this asymptotic behavior using 



the multiple shooting method 



25|. In this method one integrates the equations starting 



from the horizon towards large r, and at the same time starting from infinity towards small 
r. The two solutions should match at some intermediate point, which gives three matching 
conditions for A^, Y, U . The matching conditions should be fulfilled by adjusting the free 
parameters, but since there are only two parameters A, B in (16. ip . one cannot fulfill all three 
conditions. If one adjusts also the parameter u at the horizon, then it could be possible to 
construct global solutions, but these will exist at most for discrete sets of values of A,B,u. 
As a result, one cannot vary u continuously and so there could be no asymptotically fiat 
continuous hairy deformations of the Schwarzschild solution. 

The above argument does not exclude the U,a asymptotics (14. 7p . since they contain 
altogether 5 free parameters, which is enough to fulfill the matching conditions. 
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B. Hairy deformations of the asymptotically AdS black holes. 

Let us now choose in (15. 8 p u = C + 5u where C = C2 or C = C3 defined by (14.61) . This 
corresponds to deformations of the asymptotically AdS black holes with A(C2) = —3.0559 
or A(C3) = — 1.1812 m^. Integrating the equations shows that solutions with such boundary 
conditions exist for both values of C, provided that \5u\ is not too large. For r — )■ 00 the 
solutions approach the corresponding background black hole configuration (14. 4p . but they 
deviate from it close to the horizon. The solution profiles for C = C3 = 2.6333 are shown 
in Fig. 2, while those for C = C2 look qualitatively similar. Since the deformations do not 
change the asymptotic behavior at infinity in this case, they can be viewed as short massive 
hair localized in the horizon vicinity. 




Figure 2. Left: profiles of the asymptotically AdS solution for u = 2.8, where Nq^Qq^Yq^qq 
correspond to the background black hole (j4.4p with C = C3. Right: U'{r) for several values of u. 



It is worth noting that, unlike for asymptotically fiat solutions, a similar counting argu- 
ment does not forbid the existence of continuous deformations of the asymptotically AdS 
black holes. In order to simplify the discussion, let us set 03 = 04 = 0, in which case Eq. (l4.2p 
can be solved analytically, 

C = l±^— ^ A(C) = m2(±cosr/-l). (6.3) 
cos?7 

Let Nq^Yq^Uq be the corresponding background black hole solution (14. 4 p and consider so- 
lutions that approach this background configuration at infinity, = Ao(l + 5A), Y = 
Yq{1 + 5Y), U = Uo{l + 6U) where 6N, 6Y, 6U vanish for r 00. Inserting into (JEIl) and 
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linearizing gives 
6N 
with 

Ai = -2 



^^^ + 0i6U), 6Y = 0{6U), 



' 2 =F 5 cos rj + cos^ r] 



A, 



2 =1= 5 cos r] + cos^ r] 



(6.4) 



(6.5) 



1 =F cos 7] y 1 =F cos T] 

Since 3ft(Ai) < (unless for t] = 0, vr) and 3ft(A2) < 0, both the Ai and A2 modes are 
acceptable. The asymptotic solution therefore contains tree integration constants A, Bi,B2, 
which is enough to fulfill the three matching conditions within the shooting method. 



C. Deformations of the Schwarzschild-de Sitter black hole. 




1 




T//30 ' 




0.5 

















-0.5 
-1 


U'/2Q 


Tl/30 I 

1 1 






D 0.2 0.4 0.6 0.8 


1 1.2 

ln(r/rfe) 


1.4 1.6 1.8 



Figure 3. Deformations of the asymptotically dS solution for u = —8 (left) and u = —9 (right). 



Let us now consider the last value in ( 14. 6p . C = C4 = —8.5566, which gives rise to the SdS 
black hole (14. 4p with A = +21.5625 m^. Apart from the event horizon at r = 1, this solution 
has a cosmological horizon at r = rc where N{rc) = Q{rc) = 0. Let us set m = C + 6u. 
Integrating the equations shows that the solution becomes singular at a finite distance from 
the horizon. If 5m < then the singularity is located at r ^ rc where vanishes while Q 
does not, while the derivatives Q', [N"^)' and the energy density Tq diverge (see Fig.3). This 
implies divergence of the Riemann tensor. 

Somewhat peculiar features are shown by solutions with 6u > 0. In this case Q develops 
a simple zero at r ^ Tc but A^ remains finite, so that Tq is also finite while shows a simple 
pole. The curvature diverges at this point. Curiously, despite the latter, the solution can 
be further continued until a point where Q,N,TI^ remain finite but U',N' diverge as does 
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the energy density Tq which goes to minus infinity (see Fig.3). This produces a curvature 
singularity, as well as an infinite violation of the weak energy condition. 

Summarizing, the Schwarzschild-de Sitter black hole does not admit non-compact, regular 
hairy generalizations, since all its deformations develop a curvature singularity at a finite 
proper distance from the horizon. 

VII. GENERIC BLACK HOLES 

So far we have been considering solutions parametrically close to the background black 
holes, that is for u = Ck + Su where 6u is not too large. Let is now consider what happens for 
arbitrary u, when we deviate further and further away from the background value u = Ck- 

To begin with, solutions do not always exist, since for some values of u the argument 
of the square root in (15. 9p may become negative thus rendering the parameter ai{u) and 
the local solution (15. 8 P complex- valued. Secondly, even if ai{u) in (15. 9p is real, it should be 
positive, since we assume that A^^ grows ai r = (black hole horizon). (We do not consider 
the case where A^^ decreases at r = r/^ (cosmological horizon), however it can be treated by 
making a formal replacement N — iN , Q ^ i, Y iY, a ^ ia in the field equations.) 

Having determined the allowed values of u, we integrate the equations and find that the 
solutions always reproduce one of the types described above. They approach either one of 
the U,a backgrounds (14. 7p so that f/'(r) — )■ as r — )■ oo, or they are asymptotically AdS 
with the cosmological constant determined by (14. 6 P so that U'{r) — C2 or U'{r) C3, or 
they are compact and singular. 

We first study the upper branch solutions, with Oi = afiu) in (15. 8p . and find that they 
are non-compact only when u belongs to one of the following regions: = [—0.53; —0.65], 
I2 = [1; 1.04], = [8.3; 14.9] (see Fig.4). Solutions for u G it 1 12 are the described above 
deformations of the background AdS black hole with U'{r) = C2 and of the Schwarzschild 
solution. Solutions for u ^ are new, for r — 00 they approach the background AdS black 
hole for U'{r) = C3 but cannot be viewed as its continuous deformations, because they 
belong to the different branch and so the boundary conditions at the horizon are different. 

For all other values of u the upper branch solutions develop a curvature singularity at a 
finite proper distance from the horizon, so that they are compact. In particular, the upper 
branch solutions seem to exist for all u large and negative, but they seem to be all singular. 
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We then consider the lower branch solutions, with ai = aj~(M) in (15. 8p . and find that 
they are non-compact if u belongs to one of the four regions: If = [2.63; 3], I2 = [3.1; 29.5], 
= [29.9; 154.2], J^T = [154.4; 1876]. Solutions for u G /f are the deformations of the 
background black hole with U'{r) = C3. For u & I2 and u & we find new asymptotically 
U, a solutions, not continuously related to the Schwarzschild metric, because they belong to 
the different branch and so have different boundary conditions at the horizon. For u & 
we obtain new solutions that approach the background AdS black hole with U'{r) = C3. It 
seems that for all other values of u the solutions are compact and singular. 

When a solution from one branch is regular and non-compact, the one for the same u 
from the second branch is usually compact and singular. However, this is not a general rule, 
since, for example, the intervals and I2 have a non-zero overlap. 

The described above picture corresponds to the theory parameters m = 0.1, r/ = 1, 
C3 = 0.1, C4 = 0.3. Varying these values changes the size, position and number of intervals of 
u within which non-compact solutions exist. However, the overall picture remains the same: 
the solutions either approach asymptotically the U,a backgrounds (14. 7 P so that f/'(r) — 
as r — !■ 00, or they are asymptotically AdS with U'{r) — t- C where C is a root of (14. 6 p such 
that A(C) < 0, or they are compact and singular. If Acs + c^ > 6 then the cosmological term 
in (14. 7p changes sign and the solutions that used to be asymptotically U, a for 4c3 + C4 < 6 
become compact and singular. 

The only asymptotically flat solution we flnd for generic parameter values is the pure 
Schwarzschild black hole. Similarly, the only asymptotically dS solution is the pure SdS. 
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VIII. SPECIAL SOLUTIONS 

It is possible that new solutions could exist for special parameter values. Let us see if 
they could be asymptotically fiat. We know that the cosmological constant A(C) in ( 14. 5 p 
vanishes for C = 1. However, it will also vanish if the expression in the parenthesis in (14. 5p 
vanishes, that is for 

C = C± = — ^ f2c4 + 5c3 - 3 ± Vl2c4 + 9(C3 - l)^) . (8.1) 

Inserting this to (14. 3p . the polynomial P{C) will vanish if either C4 = — (2/3)c| or if 77 = 0, 
in which cases we obtain additional background black holes described by the Schwarzschild 
metric. Let us now see what happens if we deform these solutions by setting at the horizon 
u = C± + 6u. 

It turns out that nothing very interesting happens in the case where 77 is arbitrary and 
C4 = — (2/3)c|. There are two Schwarzschild solutions in this case, one for u = C+ and 
another one for u = C_. Choosing u = C± + 6u produces a curvature singularity at a finite 
proper distance from the horizon, similar to what is shown in Fig. 3. 

Let us now consider the case where C3, C4 are arbitrary but 77 = 0. When 77 vanishes, the 
source term in equation (12.71) vanishes too, so that the metric /^^ becomes Ricci-fiat. How- 
ever, it cannot be fiat, since our boundary conditions require that it should have a horizon, so 
that it becomes the vacuum Schwarzschild metric. Therefore, unless we change the bound- 
ary conditions, we do not recover for 7/ — the RGT theory where /^jj, is non-dynamical 
and fiat, but obtain instead the theory where f^^, is non-dynamical and Schwarzschild. The 
line element (13. ip then can be represented in the form 

dip 

where 



y(f/) = ^l-^, a{U) = AY{U) (8.3) 

with A being an integration constant. It is easy to check that this choice of Y{U),a{U) 
solves equations ( I3.6cp and ( ]3.6dp . 



Now, setting u = C± and integrating equations (15. 7p we find that the metric g^^, is also 
Schwarzschild as it should be. Setting u = C_ + 6u gives solutions which are close to the 
Schwarzschild metric in the event horizon vicinity but develop a curvature singularity at a 
finite distance from the horizon, so that this case is also not very interesting. 



21 



1.01 




01234567 012345678 

ln(r/rft) ln(r/7-h) 
Figure 5. Solutions with C3 = 0.1, C4 = 0.3, u = 1.99 and r/ = (left) and 77 = 0.02 (right). 

Qualitatively new solutions arise for u = (7+ + 6u. For example, for C3 = 0.1 and C4 = 0.3 
one has C+ = 1.97133, while setting u = 1.99 gives the result shown in Fig. 5. These solutions 
stay always close to the background Schwarzschild solution but do not tend to it for r — t- cxo 
and show instead infinitely many oscillations with a constant amplitude. One has 

N = J1-- + 5N, Q=\1-- + 5Q, U = C+x + 5U, (8.4) 
\ X \ X 

where SN, 6Q, 6U are small everywhere. Linearizing the field equations gives for large r 
SU = exp{iV2m{x + -ln{x))} + SN = -^6U+..., 6Q = - 6U + . . . , (8.5) 

2 ^/2 X 

where the real parts should be taken, the dots stand for the subleading terms, and where 
we have assumed for simplicity that 03 = 04 = 0, in which case C+ = 2. We see that the 
solutions behave as if the graviton mass was imaginary, thus providing a tachyonic version of 
the asymptotic behavior fl6.ip . fl6.2p obtained by linearizing around the C = 1 Schwarzschild 
black hole. It is interesting that the oscillations persist even for finite (but small) values of r], 
but the solutions become then asymptotically AdS and the oscillation amplitude decreases 
with r as shown in Fig. 5. 



IX. GLOBALLY REGULAR SOLUTIONS STARS AND LUMPS 

For the sake of completeness, let us briefly study what happens when there is no horizon 
and the solutions are globally regular. They can couple to a compact matter source, in 
which case they turn out to be asymptotically fiat. If the source is absent, then the solutions 
describe globally regular 'lumps of self-gravitating energy'. 
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A. Stars 

Let us return to equations fl3.6ap - fl3.6ep and restore non-zero values of p, P. We assume 
the energy density to be constant inside the star and to vanish outside, so that p(r) = if 
r < and p(r) = for r > R^. The pressure P{r) is determined by equation fl3.5p which 
should be solved simultaneously with the other equations. The pressure should vanish at 
the surface of the star, where the metrics g^^^, and /^j, should be continuous. We require 
both metrics to be regular at the origin r = 0, where the curvature should be finite. This 
leads to the local power-series solution at small r (assuming for simplicity that 03 = 04 = 0) 

N = 1+ (m2cos^r/(l --u + -u^)- — ] + 0(r^), U = ur + 0{r^), 
\ 2 2 6 / 

?/ — 1 

Y = l+m^sm'^ri x^ + Ofr^), P = p + 0(r^), Q = q + 0(r^), 

2u 

a = a{0) + 0{r^), (9.6) 

where u,p, q are free parameters which determine the value of a(0). We now wish to extend 
these local solutions numerically towards large r in order to match the fiat asymptotics 
fl6.ip .f l6l2l) . Let us again count the free parameters. We need to integrate the 5 first order 
equations - these are Eqs. fl5.7p .f l5.3cp (generalized to the case of non-zero p,P) and also 
equation (13. 5p for the pressure. In order to get the solutions within the multiple shooting 
method, we need 5 free parameters in order to match the values of 5 amplitudes A^, Y, 
U, Q, P at the matching point (the amplitude a is obtained afterwards from the algebraic 
constraint (15. 2p ). And indeed, we have in our disposal exactly 5 free parameters - these 
are u,q,p in f l9.6p and A,B in the asymptotic solution fl6.ip .f l6l2|) . We can therefore do 
the matching, and this gives us the global solutions. It is convenient to introduce mass 
functions Mg, Mf defined by iV^ ^ 1 - 2Mg/r and Y'^ = 1 - 2Mf/U, in terms of which the 
00-components of the Einstein equations (13. 2p read 

(M,)' = ^(m2cos2r/T° + p), (9.7a) 
(Mf)' = U' — m' sin" r] 7^. (9.7b) 

The typical solutions are shown in Fig.6. We see that inside the star, for r < R^, the 
pressure falls from its maximal value at the center till zero at the star surface, while the 
mass function Mg rapidly increases. The mass function Mf also increases but not as fast, 
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Figure 6. Star solutions for C3 = C4 = 0, m = 0.2, = 0.05, i?^ = 1 for three different values of 
r/. The surface of the star is located where P vanishes. The behavior of U{r) practically does not 
change when r] changes, so that the three curves for U /r almost coincide and look like as one curve. 

since it does not couple directly to p but only to 7^. Outside the star Mg decreases (for 
7] < vr/2) because it is then sourced only by Tq which is negative, while Mf still increases, 
since Tq is positive. For large r both mass functions approach the same asymptotic value 
Mg(oo) = M/(oo) = Asiv? Tj required by (16.11) . When varying 77 the parameter A almost 
does not change, so that Asiv? rj changes from the maximal value for 77 = 7r/2 to zero for 
r/ = 0. 

For 7] = tt/2 the metric g^j^ decouples and is described by the Schwarzschild solution for 

the perfect fluid. One has in this case 

^3 ^3 
Mg = p^— for r <R^ and Mg = Madm = Y ^^'^^ 

so that the mass function Mg varies only inside the star, while outside it is constant and 
equals to its asymptotic value - the ADM mass of the spacetime. The mass function Mf 
grows monotonically till the value Madm- 

If ?7 < 7r/2 then Mg grows inside the star but not as fast as for 77 = 7r/2, since the positive 
contribution of p to the right hand side of (I9.7ap is partially screened by the negative Tq . 
Outside the star Tq remains negative and continues to screen the mass of the star, so that 
Mg approaches at infinity a lesser value than it has at the surface of the star. 

For ?7 = the metric /^i, decouples and becomes flat, so that a = F = 1, Mf = 0, and 
we recover the RGT theory (unlike for the black holes). The positive mass of the star is 
then completely screened by the negative graviton energy, so that the mass function Mg 
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approaches zero at infinity. This follows from the fact that the 1/r terms in the metric 
should be absent, since the metric should now approach its asymptotic value exponentially 
fast. 



These solutions show the Vainstein mechanism of recovery of General Relativity 




Indeed, when the graviton mass m is very small, the contribution of m^Tp to the total energy 
density in Eq. fl9.7al) becomes small as compared to p. The mass function Mg then stays 
approximately constant in a large region outside the star, in which case General Relativity 
is a good approximation. 



B. Lumps of pure gravity 



Let us now return to the local solution at the origin (19. 6p and set = p = 0, in which 
case there is no source, and there is essentially only one free parameter left, u. Integrating 
towards large r, we obtain a family of solutions labeled by u and describing lumps of pure 
gravity. They have a regular center at r = 0, while for large r they behave similarly to 
the black holes and approach either the AdS backgrounds (14. 4p . or the U,a backgrounds 
(14. 7p . or develop a curvature singularity at a finite proper distance from the center. It seems 
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that the lumps are related to the black holes and can be obtained form them in the limit 
when the event horizon shrinks. Such a phenomenon is actually well known for hairy black 
holes, which can often be viewed as non-linear superpositions of a regular gravitating matter 
configuration (solution) with a vacuum black hole [26 1. 

The lump shown in the left panel of Fig. 7 approaches for r — 00 the AdS background 
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(14. 4 p with C = —0.8508, while the one on the right tends to the U, a background (14. 7p . These 
solutions can be viewed as regular deformations of the background and can be approximately 
described by replacing the constant M in (14.4p . (14.7p by a function that approaches a constant 
value at large r but vanishes at r = 0, thus insuring the regularity at the origin. 



X. CONCLUDING REMARKS 



There is something peculiar about black holes in massive gravity - all known solutions are 



of the Schwarzschild(anti)-de Sitter type 



22] ■ |13| . |14| . |15| . and it is unclear if there could 
be others, asymptotically fiat, say. New black holes could exist if only the two metrics are 
simultaneously diagonal, but since their event horizons should then coincide, this excludes 
the bimetric theories from consideration {2^. However, the argument does not apply in the 
bigravity theories, so that new black holes could exist there. Their quest was the subject of 
the above analysis. 

We did indeed find a whole zoo of new black holes with massive degrees of freedom excited. 
All of them have a regular event horizon common for both metrics, while in the outside region 
they show various types of behavior. Some of them look more conventional and support a 
short massive hair around the horizon. Some of them exhibit peculiar properties, as for 
example infinite and negative energy Tq (negative energies seem to be common for generic 
solutions), or the tachyonic features. However, none of them appear to be asymptotically 
fiat. It seems therefore that asymptotic flatness is incompatible with an event horizon as 
soon as the massive degrees of freedom are excited. This could perhaps be interpreted in 
the spirit of the no-hair theorems, many of which state that asymptotically fiat black holes 
cannot support massive hair (see {26 1 for a review and references). 

We have also constructed regular vacuum solutions without a horizon - lumps of pure 
gravity. However, they are also not asymptotically fiat and show in the far field the same 
behavior as the black holes. It seems that they can be obtained from the latter in the limit 
when the horizon shrinks to zero. 

However, adding a matter source we do find asymptotically fiat solutions which exhibit 
the Vainstein mechanism of recovery of General Relativity at finite distances. It seems 
therefore that the mechanism needs a matter source and does not work for pure vacuum 
systems like black holes. This can be supported by a qualitative argument, quite in the spirit 
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of Vainstein's original argumentation j8|. If there is a matter source, then schematically one 
has 

G{=m^T^ + T^'^'^'^, (10.1) 

and so it is clear that if the graviton mass m is very small, then one should be able in some 
way to neglect the fist term on the right as compared to the second one. This reproduces 
General Relativity. However, if T^^^^ = 0, as for the studied above black holes or lumps, 
then there is no justification for omitting m^Tj^, and General Relativity is not recovered. 
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APPENDIX. Solutions with decoupled metrics. 

In this Appendix we analyze the possibility to have a non-vanishing coefficient c in the 
metric /^i, in f l2.15p . This metric is then non-diagonal, which allows to find new solutions. 
However, these solutions are essentially of the same type as those found long ago by Isham 
and Storey 22|, in particular they reduce to those found in jl3|,jl^ when 77 — )■ 0. 

The expression r° in fl2.18p must vanish, which is possible for c 7^ if only the expression 
in the parentheses vanishes. The latter condition requires that U must be proportional to 
r, U = Cr, in which case 

rO = l{2C-3 + C3{C^ - 4C + 3) + c^{C - l)^}, (A.l) 

which can be set to zero by adjusting the value of C. Equivalently, one can consider C, C3 
as independent parameters, which will be assumed below, and this implies that 

2C-3 + C3(C^-4C + 3) 

This also implies that Tq = = X with A given by Eq. (lA.SP below, in which case the 
conservation condition (13. 3p becomes 
(9) ^,_2 _ 2(c3C - C - C3 + 2)(C^g - CQNb - aC + c'N'Q + abN) _ ^ 

(A.3) 
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Assuming for a time being that this is true, we shall solve the equations and later shall 
impose this condition on the solutions. Provided that (1A.3P is true, the energy-momentum 
tensors assume a very simple form 

Ti: = V = , (A.4) 

with 

X = iC- 1)(C3C - C - C3 + 3), A = (C3C - C3 + 2), (A.5) 

so that the field equations ( 12.6^ . ( 12171) reduce to 

= cos^ r]X6i^, (A.6) 
Ql^ = m^sm^7]~X6i; . (A.7) 

These equations describe the dynamics of the two metrics independently driven by their 
cosmological terms. The equations for g^j^^, completely decouple from those for /^^^ so that 
we can solve them independently. With gf^j^ given by fl2.17p . the solution of Eqs. (lA.6p is the 
Schwarzschild-(anti)de Sitter metric, 

q2 = = 1 - — - - cos^ . (A.8) 
r 3 

Let us now consider equations (1A.7I) for /^j,. They are slightly more difficult to solve, since 
/^jy is non-diagonal, while its components f,^^ = and f^^^ = sin^ -i? are already fixed, 
since U = Cr. However, the components /oo, for, frr are still free, because they contain 
three up to now unspecified functions a, b, c. We can consider U as the new radial coordinate, 
changing at the same time the temporal coordinate, so that t — ?■ T{t,r), r ^ U = Cr. The 
metric then becomes 

f^^dx^dx" = frr dT^ + 2fTudTdU + fuudU^ - U\d^^ + sin^ Mip^) , (A.9) 



where J'tt, fru, fuu are functions of T, U. The structure of the source term in f lA.7l) remains 



the same in the new coordinates, so that we should solve the Einstein equations with the 
cosmological term to find a metric parameterized by the radial Schwarzschild coordinate U. 
The solution is the (anti)de Sitter metric 

dip 

f^^dx^dx" = AdT^ - —--U^ (d^"^ + sin^ M^"^) , (A. 10) 
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where A{U) = 1 — ^ rri^ sin^ r] . There remains to estabhsh the correspondence between 
the T, U and t, r coordinates. Let us introduce 1-forms 

e° = VAdT, 6^ = ^, e^ = Ud^, e^ = UsmM(p, (A.ll) 

such that /^t, = TjAB^n^u ■ The correspondence between the T, U and t, r coordinates can be 



estabhshed by relating the 1-forms 9^ from (lA.lip to from fl2.13p . The two sets of 1-forms 



need not coincide but may differ by a local Lorentz rotation. This gives two conditions 

u'^ = VTT^0° + ae\ = ViT^e^ + ae^, (A.12) 

where a is the rotation parameter. These conditions explicitly read 

Cdr 



adt + cdr = Vl + a^ {y/ATdt + VAT'dr) + 



a 



Cdr 

-cN^dt + bdr = Vl + a2 + a {^/ATdt + VAT'dr). (A. 13) 



Comparing the coefficients in front of dt, dr one finds 

T{t, r)=Ct + C J f{r)dr (A.14) 



and 



/A 

Cv^TT^v^, h=VlT^^ + Ca^f, c=-Ca^ (A.15) 

VA 



with 

In order to determine the yet unspecified function a, we now use the condition that the 
expression in (lA.SD has to vanish. It will vanish if (since Q = N) 

C^N - CN% -aC + c^N^ + abN = 0, (A. 17) 

which will be satisfied if 

-A 

a = (A.18) 

2N^/A 



Together with U = Cr this finally establishes the correspondence between the t, r and T, U 
coordinates and specifies all components of f^^^. 

The above considerations give a family of the Schwarzschild-(anti)de Sitter backgrounds 
in the bigravity theory for generic values of the parameters. Some additional care should 
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be taken in order to make sure that all coefficients in the above expressions are real. For 
example, if the parameters are chosen such that / > and / < 0, then \/A will be always 
real, while will be real in the region between the black hole and cosmological horizons. 
The solutions in regions beyond the horizons can be obtained by the analytic continuation. 

When ?7 — )■ then the metric f^^ becomes fiat while the metric (7^^ does not change. The 
solutions then reduce to those obtained in the RGT theory jisl . 

Another possibility to set to zero the expression in (1A.3I) is to restrict the value of the 
coefficient C3 in such a way that the ffist factor in the numerator in (1A.3P vanishes. 



C3C-C-C3 + 2 = 0, (A.19) 



which implies that 



C-2 C'-3C + 3 J ^ , r l-C 

C3 = ^, c, = -^^-^, l = C-l, 1 = ^. (A.20) 



Since one does not need to impose the condition (1A.17P in this case, there is no equation for 



the parameter a in flA.14p . (lA.15p . flA.16|) so that it remains arbitrary. Choosing a = 0, the 
7] — limit of such solutions corresponds to the case considered in Ref . jl^ . 
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